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© Dispositif et methode pour generer une fonction ^approximation. 

© Dispositif (5) pour generer une fonction d'approximation fondee sur des premiers couples ((X,, Y,) a (Xs, 
Y 6 )) de valeurs associant une grandeur dependante <Y, a Y 6 ) a une grandeur independante (X! a Xg), et pour 
determiner des seconds couples (X A , Y' A ) de valeurs desdites grandeurs d'apres ladite fonction d'approximation. 
Le dispositif comprend : 

- des premiers moyens (10) pour determiner iterativemant au moins une fonction lineaire courante de 
regression, pour selectionner celle des fonctions lineaires courantes qui delivre I'approximation de tous les 
couples de ladite suite avec des erreurs minimales, et pour coder, a I'aide de codes specifiques (p, q), la 
fonction lineaire de regression selectionnee, 

- et des seconds moyens (17) pour determiner lesdits seconds couples (X Al Y' A ) a I'aide desdits codes 
specifiques. 

Le dispositif peut etre utilise pour le calcul de valeurs approchees de fonctions mathematiques, par exemple 
dans un reseau de neurones, ou pour determiner une fonction de regression formant une approximation de 
valeurs experimentales de mesure, par exemple des mesures dispersees issues du suivi d'un processus 
industriel. 

^invention concerne egalement une methode pour generer une fonction d'approximation. 
Application : calcul approche de fonctions mathematiques, calcul de fonctions de regression. 



r 



5 





—10 









(XA, yAi 



Fia8 



BEST AVAIUSLE COPY 



EP 0 624 847 A1 



L'invention concerne un dispositif et une methode pour generer une fonction d'approximation, celle-ci 
etant fondee sur des premiers couples de vaieurs associant une grandeur dependante a une grandeur 
independante, et pour determiner des seconds couples de vaieurs desdites grandeurs a partir de ladite 
fonction d'approximation. 

5 Un dispositif et une methode de ce type sont connus du brevet US-A- 3 789 203 qui decrit un 

generateur de fonctions operant une approximation par interpolation iterative. Ce dispositif est prevu pour 
etre utilise dans des applications de traitement de donnees necessitant un calcul de fonctions telles que 
sin(x), tg(x) par exemple. Ce dispositif ne requiert qu'une capacite de stockage minimale de la part d'un 
dispositif utilisateur. A partir de deux points appartenant a une fonction a interpoler, la methode tout d'abord 

70 interpole la fonction par une droite reliant les deux points, puis fait une approximation des ecarts entre la 
droite et la fonction par des approximations polynomials de degre croissant. Ensuite elle substitue aux 
points initiaux des points approximatifs pour reduire la longueur du segment reliant les points a traiter et 
enfin reitere les operations precedentes. 

Une telle methode necessite des ressources importantes en moyens de calcul et ne peut etre mise en 

75 oeuvre qu'avec des calculateurs performants. 

Or, il existe des applications ou une teHe methode ne peut pas etre mise en oeuvre car elles ne 
disposent pas des ressources suffisantes. De plus, pour certaines applications on peut se satisfaire d'un 
calcul approche de la fonction pour des vaieurs, en nombre limite, de la grandeur independante. 

II peut s'agir d'une fonction sigmoVde appliquee a des potentiels neuronaux delivres par au moins un 

20 neurone dans un reseau de neurones. II peut s'agir d'une autre fonction non lineaire, par exemple une 
fonction racine, pour des calculs de distances entre des etats de neurones. Les applications peuvent aussi 
concerner d'autres dispositifs comme des generateurs de fonctions, des calculateurs ou autres. 

Pour calculer une telle fonction, sans passer par une fonction d'approximation, on peut utiliser 
diffe rentes manieres. 

25 On peut effectuer le calcul mathematique exact pour chaque valeur de la grandeur independante a 
traiter, en programmant un calculateur selon les methodes connues. Une telle methode necessite d'effec- 
tuer a chaque fois les memes operations ce qui peut necessiter beaucoup de temps si le nombre de 
vaieurs est eleve. 

On peut aussi prealablement stocker dans une memoire des tables precalculees. Dans ce cas, la 
30 lecture en memoire du resultat peut etre rapide. Mais pour couvrir avec un pas assez fin toutes les vaieurs 
possibles de la grandeur independante, ii faut alors disposer de tables de grandes capacites. Ces methodes 
de calcul presentent done des inconvenients. 

D'autre part, on peut etre conduit a identifier deux grandeurs qui sont dans la dependance Tune de 
I'autre par des couples de vaieurs associant une grandeur dependante a une grandeur independante. Ainsi, 
35 dans le suivi d'un processus industriel, on peut etre conduit a mesurer par exemple un rendement R d'une 
operation en fonction de la temperature T a laquelle a ete realisee ladite opeation R = f(T). Un suivi du 
processus peut consister a consigner des lots de couples de mesures sur un graphique. Celui-ci peut etre 
utilise pour caracteriser le processus ou pour en deduire de nouveaux parametres de reglage de ladite 
operation. Ceci fait par exemple I'objet de Tarticle de H. ISHIBUCHI et H. TANAKA, "Regression analysis 
40 with interval model by neural networks" in "IEEE International Joint Conference on Neural Networks", vol. 2, 
18-21 Nov. 1991, SINGAPORE. II est souhaitable que ces nouveaux parametres soient representatifs du 
fondement de ladite operation et que les fluctuations de mesures inherentes a ce genre de processus 
soient ecartees. II est done souhaitable de determiner une approximation de la fonction f(.). 

Ainsi dans un cas il peut s'agir de mesures erratiques ou entachees d'erreurs que Ton desire 
45 representer par une fonction d'approximation. 

Dans un autre cas. on connait des vaieurs precises mais I'utilisation a en faire ne necessite pas une 
grande precision et une fonction d'approximation suffit. 

Un des buts de I'invention est de generer une fonction d'approximation avec des moyens materiels 
reduits permettant de calculer rapidement un nombre reduit de vaieurs de la grandeur dependante utiles a 
so I'application sans avoir pour cela a determiner d'autres vaieurs de la fonction d'approximation. Un but 
complementaire est de delivrer des vaieurs qui peuvent etre approchees dans la limite d'une erreur 
maximale controlee. 

Ce but est atteint avec un dispositif caracterise en ce qu'il comprend : 
- des premiers moyens : 

55 - pour determiner iterativemant au moins une fonction lineaire courante de regression en rendant 

egales, en valeur absolue, des premieres erreurs de signes alternes mesurees entre des premieres 
vaieurs de la grandeur dependante pour trois couples d'une suite desdits premiers couples, et 
respectivement des secondes vaieurs de la grandeur dependante determinees, d'apres ladite 
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fonction lineaire courante, pour les memes valeurs de la grandeur independante, 

- pour selectionner celle des fonctions lineaires courantes qui delivre I'approximation de tous les 
couples de ladite suite avec des erreurs minimales, 

- et pour coder, a I'aide de codes specifiques, la fonction lineaire de regression selectionnee, 

- et des seconds moyens pour determiner lesdits seconds couples a I'aide desdits codes specifiques. 
Ainsi avantageusement on determine une fonction lineaire de regression approchant au mieux les 

differents couples de valeurs connus. Les resultats approches ainsi • delivres forment un compromis 
satisfaisant pour de nombreuses utilisations du dispositif generateur de fonctio... 

Une fonction lineaire de regression est une fonction simplificatrice qui represente un phenomene 
complexe en reduisant les parametres significatifs. En representant la suite de couples de valeurs par des 
points dans un espace a deux dimensions, la fonction lineaire de regression devient une droite de 
regression. 

Ainsi apres avoir defini la droite de regression par des codes, on peut calculer une valeur approchee de 
la grandeur dependante en tout point de la droite de regression avec des moyens reduits pour des valeurs 
quelconques de la grandeur independante. 

L'invention concerne egalement une methode pour generer une fonction d'approximation, la methode 
comprenant : 

- une premiere phase : 

• pour determiner iterativement au moins une fonction lineaire courante de regression en rendant 
egales, en valeur absolue, des premieres erreurs, de signes alternes. mesurees entre des 
premieres valeurs de la grandeur dependante pour trois couples d'une suite desdits premiers 
couples, et respectivement des secondes valeurs de la grandeur dependante determinees, d'apres 
ladite fonction lineaire, pour les memes valeurs de la grandeur independante, 

- 1 pour selectionner celle des fonctions lineaires courantes qui delivre I'approximation de tous les 
couples de ladite suite avec des erreurs minimales, 

- et pour coder la fonction lineaire de regression selectionnee a I'aide de codes specifiques, 

- et une seconde phase pour determiner lesdits seconds couples a I'aide desdits codes specifiques. 
Les moyens mis en oeuvre par l'invention peuvent etre formes par un calculateur programme ou par un 

circuit dedie. lis peuvent aussi mettre en oeuvre des neurones. 

Un dispositif mettant en oeuvre des neurones selon l'invention peut etre utilise par un reseau de 
neurones, dont il peut notamment en constituer un sous-ensemble. En effet, pour fonctionner, le reseau de 
neurones doit disposer de moyens pour appliquer une fonction non lineaire d'activation aux potentiels de 
neurones qu'il delivre. Selon l'invention, le dispositif muni de neurones peut calculer une approximation de 
cette fonction non-lineaire d'activation. II peut egalement calculer des distances entre des etats de neurones 
en calculant une approximation d'une fonction racine carree destinee a etre exploitee dans le reseau de 
neurones. 

Lorsque la taille de la suite de couples de valeurs fournis initialement est elevee, on peut diviser la suite 
de couples en plusieurs sous-ensembles pour determiner plusieurs droites de regression et ameliorer la 
precision de I'approximation. La fonction d'approximation de la suite de couples est alors formee par une 
fonction lineaire par morceaux pour laquelle une exigeance de continuite entre les morceaux peut etre ou 
non imposee. 

Certains couples de I'ensemble de couples de valeurs peuvent avoir une influence particuliere que Ton 
peut concretiser en donnant un coefficient de ponderation specifique a chaque couple. Dans ce cas, I'erreur 
affectee a chaque couple tient compte de ce coefficient de ponderation specifique. 

Ces differents aspects de l'invention et d'autres encore seront apparents et elucides a partir des modes 
de realisation decrits ci-apres. 

L'invention sera mieux comprise a I'aide des figures suivantes donnees a titre d'exemples non limitatifs 
qui representent : 

Figure 1 : un graphique montrant une representation a deux dimensions d'un ensemble de points avec 
une droite de regression D. 

Figure 2 : un graphique montrant un ensemble de points et des droites servant a la determination d'une 
enveloppe. 

Figure 3 : un organigramme d'une premiere variante de mise en oeuvre de la methode a partir de 
triplets de points. 

Figure 4 : un organigramme d'une seconde variante de mise en oeuvre de la methode a partir de triplets 
de points. 

Figure 5 : une partie d'organigramme d'une troisieme variante de mise en oeuvre de la methode a partir 
de couples de points. 
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Figure 6 : une partie cf organigramme indiquant une preselection de points appartenant a une enveloooe 
mfeneure ou a une enveloppe superieure de I'ensemble des points enveloppe 
Figure 7 : un organigramme pour la determination des enveloppes infeneure et superieure 
F.gure 8 : un schema d'un dispositif selon Invention 

5 d?n U s r :e 9 c a s U gVn2r a ^ * ^ P0Ur " ^ * ' e ^ *"* * Session 

p'rrTa ^reltr'" diSP ° Sitif * ^ * C * H ' « * ^ ^ * Session 

Figure 1 1 : un schema d'un dispositif de transcodage pour le calcul des valeurs rt„ i a h 
>o dependante a partir d'un ensemble de droites de regression codees 9 ^ 

f2.« Un $C ^ ma d ' Une r ^ ljSa,i0n neUr ° na,e d6S moyens de comparaison de la figure 10 

! V °Z fe V* s *™™ ™ntrant une fonc.ion d'approximation formee de plusieurs drc.es 
consiuies rePreSen,a,IOnS une ^termination des limi.es de vai.dite de deux' droi.es 

>s Figure 1 5 : une representation concernant le raccordement de droites consecutives 

poin ts'TllTp" 6 p P f r ° Ximat '°" d ' une *«*on connue uniquemen, a .ravers un certain nombre de 
Z TZf n 6 ( 9Ure 1) danS Une ^sentation a deux dimensions. Chacun de ces ooints 

Pa J t P C ° UP 6 ^ Val6UrS <X ' V) re,iant ' a 9randeur -dependante x a la gn^rt^SaS v 

25 

(1) D : y = p,x + q. 

ou x et y sont des variables courantes. 

e 'iesuSe tZm^VT ^ * ^^i™ ° ^ des " reurs abs °^ s ^e er eu, 

La figure 1 represente un exemple forme de six Doints P h p Hic^il , Preaiablement. 
<s deux dimensions. A des fins d'exolication rnntl , s SP ° S8S Se ' 0n une re P r esenta.ion a 

particuiiers. ici P 3 , P A , P 5 , qu i permettent de determiner la droite de regression nntimai* „ Ui \ po,ms 
poinls p, P, ff„L p J , ' * 0 " 9 * ' e 9'»ss,on 0. qg i „ asse „i respecliwment pa, fc, s 
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Premiere variante de la premiere phase de la methode . 

La figure 3 represents la suite des etapes a mettre en oeuvre pour determiner la droite de regression 
optimale. 

5 Parmi Tensemble des points Pi... P N (bloc 100). on selectionne (bloc 102) trois points quelconques Pj, 

P jt P k avec i < j < k. Ces trois points servent a determiner la droite de regression qui minimise l'erreur sur 
y pour ces trois points. Cette determination est effectuee de maniere analytique preferentiellement par des 
moyens programmes. On determine une droite de reg.ession D pour que trois erreurs (bloc 104) : 
E PD (Pi, D). Ep 0 (P jt D), Ep 0 <P k , D) 

70 entre la droite D et chacun des points verifient : 

Ep 0 (Pi, D) = -E P0 (P h D) = Ep 0 (P k .D) 
avec 

75 

Epo (Pi. D) = Y; - (p.X, + q) 

et des relations analogues pour les autres erreurs. 

On determine la droite de regression D a I'aide des coefficients p et q de I'equation (1) tels que : 

20 



25 



30 



P = 



et 



. Y, - Xj.Yj * Xk . Y< - X,.Y k * X k . Y t - X t . Y v 

2(X k -X 1 ) 



L'erreur se rapportant a un triplet (P if P i? P k ) s'ecrit alors : 
35 E T (P h P h P k ) = \E PD (P h D)\ 

Lorsque l'erreur E T a ete ainsi calculee, on examine si les autres points de Tensemble engendrent des 
erreurs inferieures ou egales, en valeur absolue, a celles des points P*. P jt P k . Pour cela on selectionne un 
point P m (bloc 106) additionnel et on calcule la valeur absolue de l'erreur E Pm (bloc 108) apparaissant entre 

40 la valeur de la variable y au point P m et la droite D. 

Lorsque cette erreur E Pm est inferieure ou egale en valeur absolue a E T (bloc 110) (repere Y), le point 
P m additionnel est accepte et la methode se poursuit (bloc 112) avec un point additionnel suivant (bloc 
106). Si tous les points additionnels satisfont au critere | E Pm | < E T t la droite D est acceptee et ses 
coefficients sont utilises pour coder la droite de regression optimale D op , (bloc 114). 

45 Lorsque cette erreur E Pm est superieure a E T (bloc 110) (repere N), le triplet de points P*. P jf P k 
selectionne n'est pas accepte et un autre triplet de points (bloc 116) est choisi (lien 101) dans I'ensemble 
de points (bloc 102). La methode s'acheve a I'obtention de la droite O (bloc 114) satisfaisant ce critere 
meme si tous les triplets de points n'ont pas ete examines. 

II peut apparaTtre des situations ou a Tissue de I'etape 116 tous les points possibles ont ete examines 

so et aucun des triplets n'a fourni de solution (bloc 118). Dans ce cas, il est possible de reprendre le 
deroulement de la premiere variante en remplagant le test du bloc 1 10 par le test suivant : 

\E Pm \< a£ T 

55 ou a est un coefficient legerement superieur a 1. Dans cette meme situation, il est egalement possible de 
faire appel a la deuxieme variante. 

La premiere phase peut comprendre les etapes suivantes : 
A - selection de trois couples de valeurs parmi ladite suite, 
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B - calcul de la fonction lineaire courante de regression D et determination d'une erreur de triplet 

E T = \E PD \ s'y rapportant, 
C - selection d'un couple additionnel, 

D - calcul d'une erreur additionnelle E Pm entre le couple additionnel et ladite fonction, 
5 E - et lorsque | E Pm | < EK110) pour le couple additionnel. la methode reprend a I'etape C avec un 
couple additionnel suivant, 
F - et lorsque | E Pm | > E T pour au moins un couple additionnel, 

la methode reprend a Tetape A avec une selection d'un autre groupe de trois couples parmi ladite suite 
de couples, 

70 G - et lorsque \ E Pm \ < E T pour tous ies couples additionnels, la fonction lineaire courante de 
regression est codee et stockee en tant que fonction lineaire d'approximation. 
On peut choisir de scruter Tensemble des triplets en prenant un ordre croissant ou un ordre decroissant 
ou un ordre aleatoire pour effectuer cet examen. Le triplet qui sera retenu pour determiner la droite de 
regression pourra de ce fait etre detecte a un instant quelconque du deroulement de cette scrutation II 
sensuit que la rapid.te detention de la droite de regression depend de ('instant au cours duquel le triplet 
est detecte. Sa m.se en oeuvre presente un degre de complexity allant de N a N* ou N est le nombre de 
points mit.aux. Sa complexity est done reduite pour un petit nombre de points. Cette variante permet 
d obten.r une realisation materielle avec une forte parallelisation. Elle est tres peu sensible a une troncature 
des valeurs et fournit un resultat exact. 
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Deuxieme variante de la premiere phase de la methode. 

Dans cette deuxieme variante (figure 4), on selectionne successivement des triplets de points on 
calcule chaque fois une droite de regression, et, par recurrence, on selectionne celle qui deiivre rerreur'Epo 
la plus grande e'est-a-dire correspondant a I'erreur pire cas pour I'ensemble de points considers 
Selon la seconde variante, la premiere phase comprend Ies etapes suivantes : 
A - selection de trois couples de valeurs parmi ladite suite, 

B - calcul d'une fonction lineaire courante de regression D et determination d'une erreur de triplet 
t r =|£ P0 | s'y rapportant, 

C - comparaison de I'erreur E T avec une erreur optimale E op initi^see a une valeur strictement 
negative, 

D - et lorsque E T > E op , mise a jour de I'erreur optimale E op en remplagant E op par E T et mise a jour 
des codes d'une fonction lineaire optimale de regression D op en remplagant ceux-ci par Ies 
codes de la fonction lineaire courante D, 
> E - puis retour a I'etape A pour selectionner trois autres couples, 

F - et lorsque tous Ies triplets de couples de valeurs de la suite ont ete testes, Ies derniers codes de 
la fonction lineaire optimale D op constituent Ies codes de la fonction lineaire d'approximation 
: A chaque examen d'un triplet, on compare I'erreur E T dudit triplet avec I'erreur optimale precedemment 
memonsee et on met a jour I'erreur optimale E op avec la plus grande valeur de I'erreur E T determinee pour 
• chaque triplet. On met a jour egalement Ies parametres de la droite optimale. Avant la mise en oeuvre de 
I etape A, il faut initialiser la valeur E op a une valeur faible negative, par exemple -1 

Dans ce cas, la rapidite d'obtention de la droite de regression ne depend pas du mode de scrutation 
des triplets. Sa complexite est de degre N\ Une mise en oeuvre materielle peut beneficier de la grande 
regulante de I'algorithme mis en oeuvre. II est peu sensible aux troncatures des donnees et fournit une 
solution exacte. 

Troisieme variante de la premiere phase de la methode. 

Dans cette troisieme variante (figure 5), on selectionne d'abord une paire de points a laquelle on aioute 
un pom. supplemental, situe entre ces deux points, afin de former un triplet de points. Pour cela on 
mod, e les etapes A, B et C de la premiere variante. Ies autres etapes restant Ies memes. Les etapes 
modifiees sont telles que : 

Ai - modifie I'etape A en operant une selection de deux couples de valeurs appartenant a ladite 
suite, tel qu'.l existe au moins un couple additionnel intermediaire ayant une grandeur indepen- 
dante (X, - X«) comprise entre les grandeurs independantes dudit couple pour constituer au 
moins un triplet de couples, 

A2 - modifie I'etape A premierement en determinant une fonction lineaire annexe qui contient les 
deux couples selectionnes et deuxiemement en determinant des secondes erreurs entre les 
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grandeurs dependantes des couples intermediates possibles et ladite fonction lineaire annexe : 
. et lorsque ces secondes erreurs sont toutes de meme signe, selection du couple 
intermediaire fournissant la plus grande seconde erreur, en valeur absolue, pour former un 
triplet de couples de valeurs forme du couple intermediaire et des deux couples 
selectionnes, 

. et lorsque ces secondes erreurs sont de signes differents, reprise de la methode a Tetape 
A1 f 

B1 - I'etape B est effectuee avec ledit triplet selectionne, 

C1 - modifie I'etape C en selectionnant un couple additionnel dont la grandeur independante n'est 

pas comprise entre la grandeur independante des deux couples selectionnes. 
Lorsque Terreur | E Pm | est superieure a I'erreur E T (bloc 110), la methode reprend a I'etape A1 (lien 
101) avec une nouvelle selection de paire de points (bloc 102a). 

On observe que le deroulement de cette troisieme variante depend de la scrutation des valeurs et done 
des valeurs elles-memes. La complexity de la mise en oeuvre de cette variante varie entre N et N 3 ce qui 
lui confere un certain avantage par rapport aux vahantes precedentes. Une forte parallelisation des moyens 
de mise en oeuvre peut etre operee mais Implementation des moyens peut presenter un manque de 
regularity ce qui peut constituer un handicap pour une realisation integree. Cette variante est peu sensible 
aux erreurs d'arrondis des valeurs et fournit une solution exacte. 

Quatrieme variante de la premiere phase de la methode. 

Elle concerne la determination de la droite de regression a partir des enveloppes. 

II est possible de reduire le nombre de triplets a examiner en determinant des enveloppes respective- 
ment superieure et inferieure entourant les points extremes dans la representation bidimensionnelle de 
Tensemble de points. Une enveloppe superieure ou une enveloppe inferieure est definie telle qu'en joignant 
par une droite deux points adjacents quelconques de Tenveloppe, tous les autres points soient situes d'un 
meme cote respectivement de Tenveloppe superieure ou de Tenveloppe inferieure. On determine ainsi tous 
les points appartenant a ces dites enveloppes. 

La determination de la droite de regression va consister a considerer les paires de points adjacents 
d'une des enveloppes auxquels on associe un point intermediaire n'appartenant pas a ladite enveloppe 
pour constituer un triplet et operer comme cela vient d'etre decrit dans le cas des paires de points de la 
troisieme variante. Si une solution optimale n'a pas ete trouvee, on considere les paires de points de I'autre 
enveloppe. 

Pour mettre en oeuvre une enveloppe, on selectionne une paire de points adjacents appartenant a 
Tenveloppe. On determine alors s'il existe un point dispose de telle fagon que son abscisse soit 
intermediaire entre les abscisses des points selectionnes. Lorsque ce point n'existe pas on passe a une 
autre paire de points de la meme enveloppe. Pour certaines paires, lorsqu'il apparait qu'il existe un ou 
plusieurs de ces points intermediates, on choisit le point intermediaire le plus eloigne de la droite 
contenant les deux points de la paire pour former un triplet et pour determiner une droite de regression. 
Pour determiner si cette droite de regression peut etre selectionnee comme droite de regression optimale, 
la methode met en oeuvre les memes operations que celles decrites prealablement dans le cas de la 
troisieme variante. 

Pour cela on modifie la troisieme variante telle que (figure 6), prealablement a Tetape 102a (figure 5), la 
premiere phase de la methode comprend une etape (bloc 100a) de determination d'une enveloppe 
inferieure et ou d'une enveloppe superieure reunissant les points les plus extremes de Tensemble de 
points, la selection des paires de points a Tetape 102a etant faite parmi lesr points adjacents appartenant a 
Tune ou Tautre enveloppe. La selection de ladite paire de points est effectuee lorsqu'il existe au moins un 
point intermediaire ayant une abscisse situee entre les abscisses des points de la paire de points. S'il 
existe plusieurs points intermediates, on forme le triplet avec le point intermediaire le plus eloigne de la 
droite passant par les deux points qui forment la paire de points. Si une solution n'est pas trouvee avec la 
premiere enveloppe, on poursuit le traitement avec la seconde enveloppe. 

La determination des enveloppes inferieure et superieure est effectuee selon Torganigramme de la 
figure 7. L'indice des points allant croissant avec la variable x, le premier point P Q fait ainsi partie des deux 
enveloppes. Considerons d'abord Tenveloppe inferieure, les points appartenant a Tenveloppe inferieure sont 
reperes par la lettre Q. Un point courant Q v est repere par Tindice v. 

A Tetape 400 les deux premiers points de Tenveloppe sont : Q 0 = P 0 et Q, = P, . Le dernier point 
courant Q v est repere : v = 1 pour P, . Le nombre k de points traites est comptabilise. 
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s <««. to, ., ,.p, iM « f a SIS; 1X^2 ,ncre, " en ' e r,n< " ce k pour 16 w ™ p s ™ am 

droite L 0 reliant P, . P,. Le segment P, P, ~t ™£ deuxieme point Q, . On calcule la 

examine si P 3 es, au-dessus de la dro te L cl nTT,?™ ^ S69ment de '' enve,op e e bass *- On 
20 comme deuxieme segment de enJeTooe' ZZ n ^ * U Segment P * P ' est acc ^ 

examine si P, est au-dessus de d S n \S pat?e 0^""? r ^ re ' iam * « * « <* 

Pa et la droite L, sent d'abord ecartes De mJl p , T * f,9Ure *' CeCi entraTne le Point 
droite U son, ensuite auS comme ne oouva' * U d ° nC ' e P ° inf P * « la 

25 suivants (non represented). 11 P *' pu,s 1 ex 3men des points 

La complexite de mise en oeuvre de cette variant* . ^ 

depend de I'ordre de scrutation des donnees Cette comnt t T enve '° ppes est de *g* N 2 et 

Precedentesetdecefaitfournitun^ m °' ndre que Ce " e des variantes 

30 -pTss peu r bie mais 

cients de ponderation w, au poims P W n md * pendante * et "feeler des coeffi- 

faiblement variabie avec la granTeur^indan . 7 * ^ ' 0rSqUe f ° nCti ° n d ' a PP^xima«ion est 

EpoiP^D) = rV,.[V,- </>.*, + q)} 

^,r~ D - - - - >• - 

(2) p = NUMP/DET et q = NUMQ/DET 
ou les quantites NUMP. NUMQ et DET sont definies par : 

Sfl: Wi Wi(XrXi) + w i w -(x--Xi) + w,w k .(x t - Xj , 

NUMP = W,W,( Yj -Y,) + W,W,(Y k -Y,) + W,W,<Y k ' Yi ) 

Par ailleurs. I'erreur E T associee a ce triplet peut e.re exprimee e, ca.culee par : 
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avec 



55 EHP u P h P k ) = \E PD (P„D\ 

Epd{P,D) = -E PD (P h D) = E PD (P k ,D) 



ou 
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(3) E T = W ^ W y' W ^\ iX i' Y r^-yi)-(X i .Y k -X k .Y i ),(X k .Y^X i .YJl 

DET 

La premiere et la seconde variante de la methode de la premiere phase decrites precedemment 
peuvent etre mises en oeuvre en faisant intervenir les coefficients de ponderation ci-dessus. Cette mise en 
oeuvrc peut etre effectuee en programmant un calculates. 

Dans le cas ou Ton selectionne une paire de points en faisant intervenir des coefficients de ponderation 
W, la methode presente les amenagements suivants. 

Pour les points P h P,, P k qui forment le triplet, on definit un facteur F, qui lie les coordonnees 
respectives et les coefficients de ponderation respectifs des points P if P () P k . Au point central 1. on associe 
le facteur F, tel que : 



Ce facteur F, intervient dans la determination du point intermediate P, qui sera a selectionner (i < l< k) 
pour constituer le triplet en tenant compte des poids affectant les points. On choisit un point P, et on 
calcule la droite de regression Di associee au triplet P it P ( , P k et I'erreur E T , associee au triplet. 

Cette variante de la methode consiste a modifier uniquement I'etape A2 de la troisieme variante (bloc 
102b, figure 5). Cette etape determine I'existence et la valeur d'un point intermediaire servant a former un 
triplet. On cherche d'abord a former une droite de regression situee en dessous de P, et de P k Pour 
chaque point intermediaire P ( , on determine si F, = 1 et E PD (P jt D.) < 0. Si au moins un point verifie cette 
condition, il n'existe pas de droite de regression situee en dessous de Pj et P k , sinon on determine une 
grandeur G max qui est la quantite maximale parmi : 

. d'une part les quantites E PD (P jf D p ) pour tous les points intermediates, 

. d'autre part les quantites (a) suivantes. uniquement pour les points intermediates pour lesquels F, < 1 
avec 



/x F i + 1 

<«> T—[' e pd (p i' D i> • 

S'il existe au moins un point intermediaire tel que F, > 1. on definit egalement une grandeur G mjn qui 
est la valeur minimale prise par (a) uniquement pour les points intermediates tels que F, > t. On examine 
alors s'il existe au moins un point intermediaire tel que : 

Epo (Pj, D,) > 0 

G max < a E T , 

et G mio . a ^ E T , (a, coefficient > 1). S'il existe un tel point, il est choisi comme point intermediaire pour 
former le triplet. 

S'il n'existe pas de point intermediaire tel que F, > i, on teste s'il existe au moins un point 
intermediaire tel que : 
Epo (Pj, D,) £ 0 

G max S a.E T , (a : coefficient £ 1). 

Si un tel point existe, il est choisi comme point intermediaire pour former le triplet. 

Si aucun triplet n*a ete constitue, on cherche a former une droite de regression situee au-dessus des 
points Pj et P k . La meme methode est reprise en inversant le signe des erreurs E P0 . 

Si aucun point P, ne peut etre selectionne, on recommence avec une autre paire Pj, P k . 

Lorsque I'ensemble de points a traiter est trop important pour etre represents par une seule droite de 
regression, la methode code alors plusieurs droites de regression chacune etant determined selon la 
methode decrite precedemment. 

La figure 13 represente un exemple dans lequel la fonction d'approximation est constitute de plusieurs 
droites de regression. 
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Le pnncpe de la determination d'une limite optimale entre deux droits de regression adiacani^ M 

14-B pour l-exemole chST u^L l V P ' Pr6Sen,er ' 6S a " Ures presentees sur la figure 
* a pour lexemple cho.si. II sensu.t que I'approximation obtenue globalement par Dl et D2 est nntiLio 

re e ,4 Ce B e Z Ur l S ° nt V ° iSineS> C ' eSt "- dire d3nS 13 — "e«imitee par un < ail dSominT Zt 

Pour determiner la valeur d'abscisse X lim entre les deux droites ■ 
- on determine la droite Dl sur un certain nombre de points et on calcule I'erreur El maximale 
• on determ.ne la dro.te D2 sur les points restants et on calcule I'erreur E2 maximale 

- " fpSSiSS un po,nt de ,a droi,e qui — fa *■ ; - ~ - > 

■ ^ermine ,a valeur limite X lim iorsgu'i. se produi, une inversion dans le rapport entre les deux 
d'absciSef ,h ° de S ' aPP ' iqUe P " " CUrrenCe 3 P ' USieUrS dr0i,6S d * ,imi, * es ™° •«« Par P-osieurs .imi.es 

: wmMmm m 

tres faible valeur II est dm^mp h! L ° 3 9 randeur '"dependante ou < est une 

«ppllqu« pc, spbstto, ,a *L d-3 » fa ZtTrn On „S P '° CedU ' e a " ak,9l,e pet " 

**» de ,e g ,e S s«,„ tat „an s rJ,*C».l^""„Tl" * 

Dispositif cour une mise * n oeuvre de la premiere phase 
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On distingue dune part la seconde variante qui determine toutes les droites de regression passant par 
chaque combinaison realisable de triplets avec la suite de premiers couples, et les autres variantes pour 
lesquelles a chaque droite de regression determinee (associee a son erreur E T ), on examine si les points 
additionnels restants fournissent bien des erreurs additionnelles inferieures ou egales a I'erreur de triplet E T . 

La figure 9 represente un dispositif adapte a la seconde variante. II comprend : 

- une memoire 12c MEM qui stocke notamment tous les points appartenant a I'ensemble de points a 
traiter. Ces points sont represents par leurs coordonnees (x, y), et eventuellement leurs poids W ou 
leurs inverses 1/W, 

- une unite de calcul 13c COMPUT qui calcule. pour chaque triplet selectionne. la droite de regression 
adaptee a chaque triplet, c'est-a-dire les codes p. q de la droite et I'erreur E T associee au triplet 

De plus, un controleur 11 c CONTR permet d'effectuer la gestion des operations et d'adresser des 
nouveaux triplets en effectuant les lectures/ecritures de la memoire 12c et le chargement de I'unite de 
calcul 13c par de nouveaux triplets. La selection de la droite de regression qui est a conserver, c'est-a-dire 
pour cette variante celle delivrant la plus grande erreur de triplet, est effectuee par Tunite de calcul 13c 

La figure 10 correspond au cas des autres variantes pour lesquelles pour chaque droite courante de 
regression on examine si les autres points additionnels de la suite delivrent une erreur inferieure a I'erreur 
de triplet. • 

Pour cela, les premiers moyens 10 comprennent une unite de calcul 13c et une unite de comparaison 
14c formant des moyens de calcul 19c. 

L'unite de calcul 13c, transmet les codes p. q de la droite de regression du triplet courant a I'unite de 
comparaison 14c COMPAR qui determine si les points additionnels de 1'ensemble de points generent une 
erreur plus faible avec cette droite de regression que celle generee par les points du triplet courant Pour 
cela. l'unite de comparaison 14c effectue le test : 

I E Pm \> E T 

Si la droite de regression courante n'est pas acceptee (test positif). un autre triplet est selectionne et une 
procedure analogue est a nouveau effectuee. Si le test est negatif pour tous les points additionnels la 
droite de regression est acceptee et ses parametres sont charges par l'unite de calcul 13c dans la memoire 
12c. 

Pour les variantes qui operent a partir d'une paire de points, I'unite de calcul 13c est programmee pour 
determiner les points intermediates et pour former les triplets issus des paires de points, ceux-ci etant eux- 
memes eventuellement issus d'enveloppes de points. Dans ce dernier cas, l'unite de calcul 13c est 
egalement programmee pour determiner les enveloppes. L'unite de comparaison 14c effectue ensuite les 
comparaisons des erreurs se rapportant aux differents triplets. 

Pour effectuer la gestion des operations et adresser des nouveaux triplets (lorsque des triplets 
exam.nes n'ont pas donne satisfaction), le dispositif comprend, comme precedemment, un controleur 11c 
qui gere : 

- les lectures/ecritures de la memoire 1 2c, 

- le chargement de I'unite de calcul 13c par de nouveaux triplets, 

- I'appel vers I'unite de comparaison 14c de tous les points additionnels a examiner successivement 

II est possible de mettre en oeuvre la premiere phase de la methode a I'aide de moyens de 
compara.son 14c (figure 12) ayant une structure neuronale. Considerons le cas le plus complet ou il existe 
des coefficients de ponderation associes a chaque point P, Lorsque ces coefficients n'existent pas il suffit 
de leur donner une valeur unite dans les explications qui vont suivre. A litre d'exemple. on peut stocker 
dans la memoire 12c (figures 9, 10) les inverses des coefficients W ( affectes individuellement aux points P 
L unite de calcul 13c determine selon les equations (2) et (3), pour un triplet donne. les codes : -p, -q, E T et 
• Ex. 

La condition a tester dans l'unite 14c pour un point additionnel P m a tester, de parametres X m , Y m , W m , 

©st . 

Y m - pX m - (E T 'W m ) - q > 0 
ou 

Y m - p.X m + (E T /W m ) - q < 0. 
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Ces tests sont aisement realisables par une unite 14c ayant une organisation neuronale. En effet ce 
tes est construe de fonc.ions lineaires. de fonct.ons de seuillage e. de fonctions logiques aisement 
reahsees dans un reseau de neurones. 1 

Une telle unite neuronale de comparison est representee sur la figure 12. Elle comprend trois 

2 U Zl?« IT™ 65 N1 61 N2 re?0iV6nt dSS d ° nn ^ S E1 ' E2 E3 ' E4 - Le neule N3 reTo 

les sorties des neurones N1 e. N2. Chacune des entrees de ces neurones est affectee d'un coefficient 

synapt.que C, selon la technique connue mise en oeuvre dans les reseaux de neurones. Cette technique 

tZ^fsT^o^ intr0dUCti ° n t0 COmPU,in9 «" *»■ ' EEE ASSP 

Pour effectuer les tests precedents, on programme les neurones N1 et N2 selon le tableau I : 

TABLEAU I 



Neurone 


E, 


Ct 


E 2 


c 2 


E 3 


c 3 


E 4 


C* 


1 


Y m 


1 


X m 


- p 




■ Ej 


1 


-q 


2 


Y m 


1 


X m 


- p 


1A/V m 


Er 


1 


-q 



Tous les coefficients synaptiques du neurone N3 sont egaux a 1 
Les neurones N1 et N2 calculent chacun un potentiel neuronal s tel que ■ s = z- C E 

d-ac^o^est 5, H 3PPliqUer f ° nCti0n d ' aCtivation A - Pour ,es ™ et N3, la fonction 

25 A, (s) = 1 si s > 0 
Ai (s) = 0 si s < 0. 

Pour le neurone N2, la fonction d'activation A 2 est : 
30 A 2 (s) = 1 si s < 0 
A 2 (s) = 0 si s > 0. 
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La sortie du neurone N3 vaut 1 si la condition a tester est satisfaite et vaut 0 dans le cas inverse 

On observe sur le tableau I que les donnees -p. -q. -E T .+ E T sont des codes qui apparaissent comme 

%£ZJT2Z£ ' d e sch - ma de h la ,i9ure 12 Les donm§es Y - *- ,/w * so " des *™ e - S 

apparaissent sur les entrees du meme schema. Lorsque tous les points additionnels ont ete testes et 
utilises no * ** ^«™nee les codes p. q. W charges dans ,a memoire ,2 i a ors 

ut.hses pour mettre en oeuvre la seconde phase de la methode (figure 1 1 ) 

L'.nteret presente par la realisation neuronale decrite est de pouvoir paralleliser les differentes 

ZZZXZZXZ * ^ <* ^ * ~™ « «" 

secoZlZXsT S r] I 9 " 3 , 5 " 7 ** 66tMes - la seconde consistan, a calcuier des 

, 7 , ( A> VaeUfS d9S 9 randeurs est alors mise en oeuvre dans des moyens de 
decodage 17 (f.gure 11). Les codes des droites sont charges dans une memoire 12a qui au couTdela 
seconde phase est adressee, par un contro.eur Ha. pour fournir les codes des droSs de regression 
adressees. La memoire 12a. organised en .ignes par exemple. contient pour chaque die de 2 e on 

Is. S p ' q et XL °* XL est ,a lim,,e sup ' rieure d ' abscisse pour ,aa - uelle cha «- 2 ZITZ 

La memoire 12a contient ainsi un tableau de parametres 
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X L1" P C % 
correspondant aux m droites de regression stockees. 
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Les moyens de decodage 1 7 comprennent : 

- le controleur 11 a , 

- la memoire I2 a , 

- et une unite de decodage I3 a . Celle-ci regoit une valeur de requete X A de la grandeur independante 
5 pour laquelle on desire obtenir un resultat Y' A issu de la fonction d'approximation. Pour cela le 

controleur 1 1 a presente successivement certaines lignes de codes de la memoire 12 a et, pour chaque 
ligne de codes, I'unite de decodage 13 a teste le code x L pour determiner si la valeur X A de la 
grandeur d'entree peut etre traitee par la droite de regression codee par ladite ligne. Si cela n'est pas 
le cas, I'unite 13 a requiert du controleur 11 a une autre ligne de codes. Lorsque cette condition est 
w verifiee, I'unite 13 a effectue le calcul : 

Y' A = p n .X A + q n 

si la ligne d'ordre n est la ligne adequate. 
75 Comme cela vient d'etre decrit, les premiers moyens 10 peuvent constituer un dispositif de codage et 
les seconds moyens 17 peuvent constituer un dispositif de decodage, les deux dispositif etant separes. lis 
peuvent en effet etre distants Tun de I'autre, les operations de codage etant separees des operations de 
decodage. Mais les deux moyens 10 et 17 peuvent aussi former un dispositif de codage/decodage unique. 
Dans ce cas certains organes peuvent etre affectes successivement a plusieurs operations. Notamment : 
20 - il peut s'agir d'une memoire 12 unique qui peut stocker les parametres des couples de valeurs et les 
codes des droites de regression determinees. Les unites 12a et 12c forment alors la memoire 12, 

- il peut s'agir d'un controleur 1 1 unique qui regroupe les fonctions des controleurs 1 1c et 11a, 

- I'unite de calcul 13c et I'unite de decodage 13a peuvent egalement former une seule unite 13. 

It est egalement possible de realiser les moyens de decodage 17 (figure 11) selon une organisation a 
structure parallele. Une solution consiste alors a comparer, en parallele, la grandeur de requete X A a tous 
les codes x L . On obtient un ensemble de signaux T L representant les resultats de ces tests de comparaison. 
On realise ensuite, en parallele, pour toutes les valeurs de L, un test qui combine deux a deux les signaux 
T L et T L+t et qui indique si la grandeur de requete X A est dans I'intervalle d'indice L Ce test est verifie pour 
une unique valeur de L, notee n. On accede alors aux codes p n et q n correspondants, et on en deduit la 
grandeur Y' A de la meme maniere que precedemment. 

Les premiers moyens 10 (de codage) associes aux seconds moyens 17 (de decodage) peuvent etre 
utilises pour determiner la valeur d'une fonction d'approximation par au moins une droite de regression. 
Ces determinations peuvent etre effectuees pour des valeurs quelconques de la grandeur independante (a 
I'mterieur de bornes predefines definissant I'etendue de Taction de chaque droite de regression). Cette 
methode evite de stocker inutilement des tables de valeurs, pour lesquelles toutes les valeurs ne seront pas 
exploiters. Selon la methode on ne determine que les seules valeurs necessaires a Implication. La 
methode selon Invention presente I'interet de ne calculer que les valeurs necessaires. Le dispositif de 
{'invention peut etre utilise en combinaison avec un processeur neuronal pour calculer une approximation 
d'une fonction non lineaire par exemple une fonction sigmoide. Un processeur neuronal calcule des 
potentiels neuronaux issus de sommes de produits de coefficients synaptiques par des etats de neurones. 
La technique est developpee dans le document de R.P. LIPPMANN deja cite. Chaque potentiel neuronal 
doit etre soumis a Taction d'une fonction non lineaire. Pour cela, chaque potentiel neuronal constitue la 
grandeur independante decrite precedemment et le calcul de la fonction th(x) par exemple est effectue pour 
les seules valeurs de requetes utiles. La methode fournit des valeurs approximees. mais cette approxima- 
tion peut etre obtenue avec une precision qui peut etre accrue selon le nombre de droites de regression 
calculees et eventuellement les coefficients de ponderation W. La precision reste neanmoins tributaire de la 
precision des coordonnees des couples de depart. 

Une telle methode est particulierement interessante pour les calculs de fonctions connues (telles que 
des fonctions mathematiques) ou de fonctions inconnues explicitement, par exemple une fonction represen- 
tee par des points de mesures, que Ton desire simplifier par une fonction lineaire de regression. 

Dans sa version neuronale, I'invention est interessante dans des applications neuronals car elle 
apporte non seulement un traitement homogene mais aussi une grande compacite de ('architecture 
materielle necessaire. 

II est possible de combiner les variantes de la methode. dans leurs mises en oeuvre materielles pour 
combiner leurs avantages. Ainsi par exemple. les combinaisons suivantes sont avantageuses : 
Dans le cas ou les poids existent, deux combinaisons apparaissent interessantes : 
. Dans la premiere combinaison, on utilise tout d'abord la premiere variante basee sur des triplets de 
points. Si une solution n'est pas trouvee. on enchaine avec la deuxieme variante basee aussi sur des 
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triplets de points. Ceci permet d'obtenir les avanlages suivants : 

Quelles que soient les donnees. on obtient une solution. 
. Dans le cas ou on exige qu'on obtienne toujours une solution exacte. cette combinaison permet 
d'obtenir la rap.d.te maximale. En effet. dans le cas general la solution est fournie directement par 
la premiere vanante appelee, qui est celle fournissant le plus rapidement une solution exacte dans 
le cas ou les poids existent. La deuxieme variante est beaucoup plus lente mais elle permet de 
garantir que la combinaison fournit toujours une solution. 
. Dans la deuxieme combinaiscn. on utilise tout d'abord la variante basee sur des paires de points Si 
elle ne trouve pas une solution (qui peut etre approchee), on utilise la deuxieme variante basee sur 
des triplets de points. Ceci permet d'obtenir les avantages suivants : 
. Quelles que soient les donnees. on obtient une solution (eventuellement approchee). 
. Cette combinaison est en moyenne plus rapide que la precedente. 
. En contrepartie. elle ne fournit generalement qu'une solution approchee. 
Dans le cas ou les poids n'existent pas. une combinaison interessante consiste a utiliser tout d'abord la 
variante basee sur I'enyeloppe des points. Si elle ne trouve pas de solution, on utilise la deuxieme variante 
basee sur des triplets de points. Ceci permet d'obtenir les avantages suivants : 
. Quelles que soient les donnees, on obtient une solution, 
. Cette combinaison permet d'obtenir la meilleure rapidite moyenne 

Lorsque la fonction lineaire de regression a ete determinee selon la methode decrite, il est possible de 
determiner non seulement des grandeurs Y' A ressortissan. a ladite fonction mais egalement des grandeurs 
ressortissan a la denvee de ladite fonction. En effet. pour une valeur de requete X A de la grandeur 
independante. ,l suffit d'extraire des moyens de memorisation la valeur p ressortissant a la valeur X A pour 
connaitre la valeur de cette derivee correspondante. 

Revendications 

1. Dispositif (5) pour generer une fonction d'approximation fondee sur des premiers couples ((X, Y,) a 
(Xg Y 6) ) de valeurs associant une grandeur dependante (Y, a Y 6 ) a une grandeur independante (X, a 
Xs), et pour determiner des seconds couples (X A . Y' A ) de valeurs desdites grandeurs d'apres ladite 
fonction d approximation caracterise en ce que le dispositif comprend : 

- des premiers moyens (10) : 

- pour determiner iterativemant au moins une fonction lineaire courante de regression en 
rendant egales. en valeur absolue. des premieres erreurs, de signes alternes. mesurees entre 
des premieres valeurs (Y 3 , Y.. Y 5 ) de la grandeur dependante pour trois couples <X 3 Y 3 ) (X* 
1. { ZV J, ] d ' UnS SU ' te d6SditS P remiers tuples, et respectivement des secondes valeurs 
(Y 3 , Y«. Ys) de la grandeur dependante determines, d'apres ladite fonction lineaire courante 
pour les memes valeurs (X 3 , X» . X 5 ) de la grandeur independante 

- pour selectionner celle des fonctions lineaires courantes qui delivre ('approximation de tous les 
couples de ladite suite avec des erreurs minimales, 

- et pour coder, a I'aide de codes specifiques '(p. q), la fonction lineaire de regression 
selectionnee. 

cod d eTspSiqu e r yenS (17> P ° Ur d4 ' erminer ,eSdi,s seconds couples Y '*> ■ rai de desdits 

2. Dispositif selon la revendication 1 caracterise en ce que les premiers moyens (10) comprennent ■ 

- des moyens (12c) pour memoriser les grandeurs des premiers couples et lesdits codes specifi- 
ques, ^ 

- des moyens (13c) (19c) pour calculer successivement les codes specifiques (p, q) de fonctions 
lineaires courantes de regression et leurs premieres erreurs, egalisees en valeur absolue 
respectives (E T ) et pour selectionner au moins une fonction lineaire courante de regression 

" calculer" 10 ^"* ° ^ C °" tr6ler ' 6S moyens (12c) pour ^moriser et les moyens (13c) pour 

3. Dispositif (5) selon les revendications 1 ou 2 caracterise en ce que pour selectionner ladite fonction 
lineaire, les premiers moyens comprennent des moyens (14c) : 

- PO^ determiner pareillement des erreurs additionnelles pour les autres couples (X, Y,) (X 2 Y ? > 
(Xs , Y 6 ) de valeurs de ladite suite. ' 

- pour comparer les erreurs additionnelles avec lesdites premieres erreurs egalisees (Et), 
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- et pour seiectionner la fonction lineaire courante qui fournit des erreurs additionnelles inferieures 
ou egales en valeur absoiue auxdites premieres erreurs egalisees (E T ). 

Dispositif selon la revendication 2 caracterise' en ce que, pour seiectionner la fonction lineaire de 
regression, les moyens pour calculer (13c) selectionnent celle des fonctions lineaires courantes de 
regression qui delivre des premieres erreurs egalisees (E T ) maximales. 

Dispositif selon une des revendications 1 a 4 caracterise en ce que les seconds moyens (17) 
comprennent : 

- des moyens (12a) pour memoriser lesdits codes specifiques d'au moins une fonction lineaire de 
regression, 

- des moyens (13a) pour seiectionner, sur requete d'une valeur de grandeur independante (X A ). 
une fonction lineaire de regression ressortissant a ladite valeur, et pour calculer, a partir des 
codes specifiques regus des moyens pour memoriser, une valeur de grandeur dependante (Y' A ) 
associee a la valeur de grandeur independante (X A ), 

- et des moyens (1 la) pour adresser aux moyens pour calculer, et ur leur demande, les codes 
specifiques des fonctions lineaires de regression. 

Dans un dispositif pour generer une fonction d'approximation fondee sur des premiers couples ((X, , Yi ) 
- (Xg, Y 6 )) de valeurs associant une grandeur dependante (Yi - Y 6 ) a une grandeur independante (Xi - 
Xg), methode pour generer ladite fonction d'approximation a partir des premiers couples (X; . Y,) (Xe, 
Y 6 ) de valeurs, et pour determiner des seconds couples (X Al Y' A ) de valeurs desdites grandeurs 
d'apres ladite fonction d'approximation, caracterisee en ce que, la methode comprend : 
,- une premiere phase : 

- pour determiner iterativement au moins une fonction lineaire courante de regression en 
rendant egales, en valeur absoiue, des premieres erreurs E PD , de signes alternes, mesurees 
entre des premieres valeurs (Y 3l Y 4 , Y 5 ) de la grandeur dependante pour trois couples (X 3 , 
Y 3 ), (X4, Y 4 ) (X 5 , Y 5 ) d'une suite desdits premiers couples, et respectivement des secondes 
valeurs (Y* 3 , Y f 4, Y' s ) de la grandeur dependante determined, d'apres ladite fonction lineaire, 
pour les memes valeurs (X 3 , X4, X 5 ) de la grandeur independante, 

- pour seiectionner celle des fonctions lineaires courantes qui delivre I'approximation de tous les 
couples de ladite suite avec des erreurs minimales, 

- et pour coder la fonction lineaire de regression selectionnee a I'aide de codes specifiques. 

- et une seconde phase pour determiner lesdits seconds couples (X A , Y' A ) a I'aide desdits codes 
specifiques. 

Methode selon la revendication 6 caracterisee en ce que la premiere phase comprend les etapes 
suivantes : 

A - selection (102) de trois couples (P3, P4. P5) de valeurs parmi ladite suite, 

B - calcul (104) de la fonction lineaire courante de regression D et determination d'une erreur de 

triplet E r =|fpo| s'y rapportant, 
C - selection (106) d'un couple additionnel (P1), (P2), (P6), 

D - calcul (108) d'une erreur additionnelle E Pm entre le couple additionnel et ladite fonction, 

et lorsque | E Pm \ < EK110) pour le couple additionnel, la methode reprend a I'etape C avec 
un couple additionnel suivant. 

et lorsque | E Pm | > E T pour au moins un couple additionnel (110), la methode reprend a 
I'etape A avec une selection d'un autre groupe de trois couples parmi ladite suite de 
couples, 

G - et lorsque \ E Pm \ < E T pour tous les couples additionnels (112), la fonction lineaire courante 
de regression est codee et stockee (1 14) en tant que fonction lineaire d'approximation. 

Methode selon la revendication 6 caracterisee en ce que la premiere phase comprend les etapes 
suivantes : 

A - selection (102) de trois couples (P3, P4. P5) de valeurs parmi ladite suite. 

B - calcul (104) d'une fonction lineaire courante de regression D et determination d'une erreur de 

triplet E T =\E PD \ s'y rapportant, 
C- comparaison (210) de I'erreur E T avec une erreur optimale E op initialised a une valeur 

strictement negative. 
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et lorsque E T > E op , mise a jour (1 14) de I'erreur optimale E op en remplacant E op par E T et 
mise a jour des codes d'une fonction lineaire optimale de regression D op en rempiagant 
ceux-ci par les codes de la fonction lineaire courante (D), 
puis retour (116) a I'etape A pour selectionner trois autres couples (102), 
et lorsque tous les triplets de couples de valeurs de la suite ont ete testes (116), les derniers 
codes de la fonction lineaire optimale D op constituent les codes de la fonction lineaire 
d'approximation (114). 

9. Methode selon la revindication 8 caracterisee en ce que la premiere phase comprend les etapes 
modifiees suivantes : 

A1 - modifie I'etape A en operant une selection (102a) de deux couples de valeurs appartenant a 
ladite suite, tel qu'il existe au moins un couple additionnel intermediate ayant une grandeur 
independante (Xi - Xg) comprise entre les grandeurs 
independantes dudit couple pour constituer au moins un triplet de couples de valeurs. 

A2- modifie I'etape A premierement en determinant (102b) une fonction lineaire annexe qui 
contient les deux couples selectionnes et deuxiemement en determinant des secondes 
erreurs entre les grandeurs dependantes des couples intermediates possibles et ladite 
fonction lineaire annexe : 

. et lorsque ces secondes erreurs sont toutes de meme signe, selection du couple 
intermediate fournissant la plus grande seconde erreur, en valeur absolue, pour 
former un triplet de couples de valeurs forme du couple intermediate et des deux 
couples selectionnes, 

. et lorsque ces secondes erreurs sont de signes differents, reprise de la methode a 
i'etape Al, 

B1 - I'etape B (104) est effectuee avec ledit triplet selectionne, 

C1 - modifie I'etape C (106) en selectionnant (106a) un couple additionnel dont la grandeur 
mdependante n'est pas comprise entre la grandeur independante des deux couples selec- 
tionnes. 

10. Methode selon la revendication 9 modifiee en ce que. prealablement a I'etape A1, lesdits couples de 
valeurs etant represents par des points (Pi - P6) dans un espace a deux dimensions, la premiere 
phase de la methode comprend une etape (100a) de determination d'une enveloppe inferieure et/ou 
d'une enveloppe superieure reunissant les points les plus extremes de ('ensemble de points la 
selection des deux couples de valeurs a I'etape A1 etant faite parmi les points adjacents appartenant a 
I une ou I'autre enveloppe. 

11. Methode selon une des revendications 6 a 10 caracterisee en ce qu'a chaque valeur, prise dans la 
suite de valeurs de la grandeur independante, est associee un coefficient specifique de ponderation qui 
pondere lesdites premieres erreurs. 

12. Methode selon une des revendications 6 a 11 caracterisee en ce qu'elle genere une fonction 
d approx.mation comprenant une suite de fonctions lineaires de regression raccordees entre elles. 
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